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1. Introduction 

Let (Xj, Yj), 1 < z < be a sample of size iV > 1 from a couple (X, Y) of in¬ 
dependent random variables (rv’s) defined over some probability space (fl, A, P), 
with continuous marginal distribution functions (df’s) F and G respectively. Suppose 
that X is truncated to the right by Y, in the sense that Xj is only observed when 
Xi < Yj. This model of randomly truncated data commonly hnds its applications 
in such areas like astronomy, economics, medicine and insurance. In the actuarial 
world, for instance, it is usual that the insurer claim data do not correspond to 
the underlying losses, because they are truncated from above. Indeed, when fac¬ 
ing large claims, the insurance company specifies an upper limit to the amounts 
to be paid out. The excesses over this hxed threshold are then covered by a rein¬ 
surance company. This kind of reinsurance is called excess-loss reinsurance (see, 
e.g., Embrechts et ah, 1997). In life insurance, the upper limit, which may be ran¬ 
dom, is called the cedent company retention level whereas in non-life insurance, it 
is called the deductible. The usefulness of the statistical analysis under random 
truncation is shown in Herbst (1999), where the author applies truncated model 
techniques to estimate loss reserves for IBNR (incurred but not reported) claim 
amounts. For a recent discussion on randomly right-truncated insurance claims, one 
refers to Escudero and Ortega (2008). Some examples of truncated data from as¬ 
tronomy and economics can be found in Woodroofe (1985) and for applications in 
the analysis of AIDS data, see Wang (1989). In reliability, a real dataset, consisting 
in lifetimes of automobile brake pads and already considered by Lawless (2002) in 
page 69, was recently analyzed in ? as an application of randomly truncated heavy¬ 
tailed models. Since the focus is on datasets that contain extreme values, then it 
would be natural to assume that both survival functions F := 1 — F and G := 1 — G 
are regularly varying at inhnity with tail indices 71 > 0 and 72 > 0 respectively. 
That is, we have, for any a; > 0, 

Ito and lim SM = (i.i) 

2 ->-oo F G {z) 

This class of distributions, which includes models such as Pareto, Burr, Frechet, 
Levy-stable and log-gamma, takes a prominent role in extreme value theory. Also 
known as heavy-tailed, Pareto-type or Pareto-like distributions, they provide appro¬ 
priate descriptions for large insurance claims, log-returns, large price fluctuations, 
etc... (see, e.g.. Resnick, 2006). 
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Let us denote (Xj,Fj), i = 1,n to be the observed data, as copies of a couple 
of rv’s (X, Y) , corresponding to the truncated sample (Xj, Yj) , i = 1,N, where 
n = uat is a sequence of discrete rv’s. By of the law of large numbers, we have 
n/N —)■ p := P (X < Y) as n —)■ oo. For convenience, throughout the paper, the 
convergence in probability of n and/or any of its subsequences is simply denoted by 

p 

instead of —)■ . The joint distribution of Xi and Y* is 
H{x,y) ■.= P{X <x,Y <y) 

= P (X < X, Y < p I X < Y) = Tf (min (x, z)) dG (z). 

Jo 

The marginal df’s of the observed X's and Y's, respectively denoted by F and G, 
are equal to 

F (x) := p-^ [ G (z) dF (z) and G (y) := p'^ T F (z) dG {z) . 

Jo Jo 

It follows that the corresponding tails are 

poo poo 

F (x) = -p-^ / G (z) dF (z) and G (y) = / F (z) dG (z). (1.2) 

J X J y 

It is clear that the asymptotic behavior of F simultaneously depends on G and F, 
while that of G only relies on G. Making use of Proposition B.1.10 in de Haan and Ferreira 
(2006), for the regularly varying functions F and G, we may readily show that both 
G and F are regularly varying at infinity as well, with respective tail indices 72 and 
7 := 'yi'^ 2 / (71 + 72 ) • If is worth noting that the issue of analyzing extreme values 
in the context of random truncation, is at an early stage. Indeed, the first contribu¬ 
tion was made in the recent paper of Gardes and Stupfler (2015), where the authors 
exploited the above relation between the three indices to dehne an estimator of 71 
by considering the classical Hill estimators of 7 and 72 (Hill, 1975) as functions of 
two distinct numbers of top statistics. However, they did not handle the case where 
these numbers are equal because of the difficulty in assessing the dependence be¬ 
tween the two Hill estimators. In the present work, we introduce a tail product-limit 
process for which we provide a weighted Gaussian approximation as well. This tool 
will be very helpful when dealing with the estimation of any tail related quantity. 

In particular, it will lead to the asymptotic normality of the extreme value index 
estimator that we define, under random right-truncation, as a function of a single 
sample fraction of upper order statistics. But, prior to describing our estimation 
methodology, let us note that, as mentioned by Gardes and Stupfler (2015), in or¬ 
der to ensure that it remains enough extreme data for the inference to be accurate. 
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we need to impose the condition 71 < 72 - In other words, we consider the situation 
where the tail of the rv of interest X is not too contaminated by the truncation rv 
Y. Since F and G are heavy-tailed, then their right endpoints are infinite and thus 
they are equal. Hence, from Woodroofe (1985), we may write 

L F(2/) X C{y)' 

where 

C{z) :=P{X < z <Y) = F{z)-G{z) . (1.4) 

Differentiating (1.3) leads to the following crucial equation 

C (x) dF (x) = F (x) dF (x), (1.5) 


(see, for instance, Strzalkowska-Kominiak and Stute, 2009), whose solution is de¬ 
fined by F (x) = exp — dF (z) /C (z). Replacing F and C by their respective 
empirical counterparts 


Fn (x) := n ^ 1 (Xj < x) and C„ 


[x := n 


^ 1 (X, < X < Yi ), 


Z=1 


i=l 


yields the product-limit estimator 


F„ (i) := TT exp I — ( \ , 


to the underlying df F. The hrst mathematical investigation on this estimator may 
be attributed to Woodroofe (1985) and the central limit theorem under random 
truncation was established by Stute and Wang (2008). Note that the approximation 
exp (—t) 1 — t, for small f > 0 , results in the well-known estimator introduced by 
Lynden-Bell (1971). Let us now introduce a tail product-limit process corresponding 
to F„ as follows: 


i>0. (1.6) 

\ -t* n \^n—k:n) / 

where Xi:„ < ... < Xi:„ denote the order statistics pertaining to Xi,...,X„ and 
/c = is a sequence of discrete rv’s satisfying 


l<A;<n, k ^ 00 and k/n —)■ 0 as n —)■ cx). (1.7) 

Observe that, in the case of complete data we have F„ = Fn with Fn (Xn-k-.n) = k/n 
and thus the process defined in ( 1 . 6 ) becomes 

Dn (x) := \/k (^^Fn (xXn-k-.n) “ X~^Fl'^ _ 
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By jointly applying Theorems 2.4.8 and 5.1.4 (pages 52 and 161) in de Haan and Ferreira 
(2006) we have that, for a:o > 0 and 0 < ^ < 1 / 2 , 


sup X 

X>X0 


( 1 / 2 - 0/71 


Dn {x) — r (x; W) — X 


-1/71 


X 


'rl/71 


Tl'Jl 


-VkA(j (n/k) 


4o, 


( 1 . 8 ) 


provided that F fulfills the second-order regular variation condition with auxiliary 
function Aq tending to zero, not changing sign near infinity, having a regularly 
varying absolute value with index ti < 0 and satisfying ^/HAq {n/k) = O (1). Here 
r (x; W) := W (1) with {W (s); 0 < s < 1} being a standard 

Wiener process. Many authors used this approximation to establish the limit dis¬ 
tributions of several statistics of heavy-tailed distributions, such as tail index esti¬ 
mators (see, e.g., de Haan and Ferreira, 2006, page 76) and goodness-of-£t statistics 
(Koning and Peng, 2008). The main goal of this paper is to provide an analogous 
result to ( 1 . 8 ) in the random truncation setting through the tail product-limit pro¬ 
cess ( 1 . 6 ) , which, to the best of our knowledge, was not addressed yet in the extreme 
value theory literature. 


The rest of the paper is organized as follows. In Section 2, we present our main result 
which consists in a Gaussian approximation to the tail product-limit process D„ (x). 
As an application, we introduce, in Section 3, a new Hill-type estimator (Hill, 1975) 
for the tail index 71 and we establish its consistency and asymptotic normality. The 
hnite sample behavior of the proposed estimator is checked by simulation in Section 
4. The proofs are postponed to Section 5 and some results that are instrumental to 
our needs are gathered in two lemmas in the Appendix. 


2. Main results 


Weak approximations of extreme value theory based statistics are achieved in the 
second-order framework (see de Haan and Stadtmliller, 1996). Thus, it seems quite 
natural to suppose that df’s F and G satisfy the well-known second-order condition 
of regular variation that we express in terms of the tail quantile functions. That is, 
we assume that for x > 0 , we have 


and 


lim 

t^OO 


lim 

t^OO 


Up {tx) /Up (f) — x'^i 
Ap (f) 

Ug (tx) /Ug (t) - x^2 


= x'^^ 


X 


Tl 


= x'*'^ 


X 


(2.9) 


( 2 . 10 ) 


Ag ( 1 ) T2 

where ri, r 2 < 0 are the second-order parameters and Ap, Ag are functions tending 

to zero and not changing signs near infinity with regularly varying absolute values 
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at infinity with indices ri, T 2 respectively. For any df K, the function (t) := 
(1 — stands for the tail quantile function. 


Theorem 2.1. Assume that both second-order conditions (2.9) and (2.10) hold with 
7i < 72. Let k = kn be a sequence satisfying (1.7), then there exist a function 
Aq (t) ~ Af (1/F (Uf (f))) and a standard Wiener process {W (s); 0 < s < 1} , 
defined on the probability space ( 12 , A, P) , such that, for 0 < f < 1/2 — 7/72 and 
Xq > 0 , we have 


sup 

X>XQ 


(a:)-r(a;;W) 


— X 


•rihl _ 

- VkAo in/k) 

llTl 


4o, 


as n —)■ 00 , provided that \/kAo{n/k) = 0(1), where {r(x;W); a; > 0} is a 
Gaussian process defined by 


r(x;W) ;= —{x^/^W(x-^/^) - W(l)} 
7i 

/•i 


+ 


7 


/ S-O 72-1 (x-^/^s) - W {s)}ds. 


7i + 72 

Remark 2.1. A very large value of 72 yields a '-y-value that is very close to 71 , 
meaning that the really observed sample is almost the whole dataset. In other words, 
the complete data case corresponds to the situation when I /72 = 0 , m which case we 
have 7 = 7 i. It follows that 

r ^ 

_1 — / s- 7 / 72 - 11 ^ 1 / 7 ^ (x-iOs) - W (s) I ds = 0, 

71 + 72 70 ^ ^ ^ ^ 

and therefore T (x; W) = W (x“^0^) — x“^O^W (1), which agrees with the weak 

approximation ( 1 . 8 ). 


3. Tail index estimation 

We start the construction of our estimator by noting that from Theorem 1.2.2 in 
de Haan and Ferreira (2006), the hrst-order condition (1.1) (for F) implies that 

1 /-oo _ 

lim =— / x“^F (x) dx = 71 , 


which, by an integration by parts, becomes 


1 f°° X 

lim =— / log—dF(x) = 71 . 
i^ooF(f)A ^t ^ ^ 


(3,11) 


Replacing F by F„ and letting t = Xn-k-n yields 


7i := = 


Fn iyXji_k-,f 


log 


X 


'X, 


n — k\n 




n—k:n 


-dYn (a:), 
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as a new estimator to 71 . By setting (x) := 1 {x > Xn-k-.n} ^og{x/Xn-k:n) and 
(x) := 1 {x > Xn-k-.n} , this may be rewritten into 

/“ (x) dFn (x) 


7i = 


20 _ 

(x) dFn ( 2 :) 


From the empirical counterpart of equation (1.5) we get 

{Xn 


poo 1 _ _ 

/ (x) dFn {x) = - 

^ ^^-k 


Cn {Xnr 


and 


poo -1 

/ ^nHx)dFn{.x) = - 
do " 

Finally, changing z to n — z + 1 yields 


log {^Xnn/Xn—k-.n) ) 


Fn {Xi,n) 


i=n—k 


CniXnnY 


7i = 


E 


Ffi (^Xn—i-{-l:n) 


-1 


E F^ yXn—i-\-l-,n) 1 ^n—i+l:n 

TrTTT-7 log 


Cn (^n—i+l:n) j Cn (^n—i+l:n) Xn—k.n 


Remark 3.1. For complete data, we have Fn=Fn=Cn and consequently reduces 
to the classical Hill estimator (Hill, 1975). 


Theorem 3 . 1 . Assume that (1.1) holds with 71 < 72 and let k = kn he an integer 
sequence satisfying (1.7). Then 71 —)■ 71 in probability. Assume further that both 
second-order conditions (2.9) and (2.10) hold and VkAo [n/k) = O (1) , then 

V 7 ( 7.-70 = #^- 7 W( 1 ) 

1 - Ti 

H- - — [ (72-71-7 log s) ( 2 s + Op (1). 

7 i + 72 Jo 

Corollary 3 . 1 . If, in addition to the assumptions of Theorem 3 . 1 , we suppose that 
VkAo {n/k) A, then 

—, cr^ ) , as n ^ 00, 

vi-^i / 

where 

:= V (1 + 71/72) (1 + (71/72)^) / (1 - 71/72)^ • 

Remark 3 . 2 . In the case of complete data we have, from Remark 2.1, cr^ = 7^. 
It follows that Vk (71 — 71) ^ A/" (A/ (1 — ri), 7^), as n —)■ 00, which meets the 
asymptotic normality of the classical Hill estimator (Hill, 1975 ), see for instance. 
Theorem 3 . 2.5 in de Haan and Ferreira ( 2006 ). 
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4. Simulation study 

This study, just intended for illustrating the performance of our estimator, is realized 
through two sets of truncated and truncation data, both drawn from Burr’s model: 

F (x) = (l + , G (x) = (l + x^/*^) , a: > 0, 

where h, 71,72 > 0. The corresponding percentage of observed data is equal to 
V = 72/(71 + 72 )- We £x 5 = 1/4 and choose the values 0.6 and 0.8 for 71 and 70%, 
80% and 90% for p. For each couple ( 71 , p), we solve the equation p = 72/(71 + 72 ) 
to get the pertaining 72 -value. We vary the common size N of both samples 
(Xi,...,X 7 v) and (Yi,..., Y^r), then for each size, we generate 1000 independent 
replicates. Our overall results are taken as the empirical means of the results ob¬ 
tained through all repetitions. To determine the optimal number (that we denote 
by k*) of upper order statistics used in the computation of 71 , we apply the algo¬ 
rithm of (Reiss and Thomas, 2007, page 137). The performance of the newly dehned 
estimator, in terms of absolute bias and root of the mean squared error (rmse) is 
summarized in Table 4.1, where we see that, as expected, the size of the initial 
sample influences the estimation: the larger iV, the better the estimation. On the 
other hand, we note that the estimation accuracy decreases when the truncation 
percentage increases, which seems logical. Finally, we observe that the estimation 
of the larger value of the tail index is less precise. 

5. Proofs 

5.1. Proof of Theorem 2.1. Set Ui := F {Xi) and dehne the corresponding uni¬ 
form tail empirical process by (s) := Vk (Un (s) — s), for 0 < s < 1 , where 
(s) := k~^ < ks/n). The weighted weak approximation to (s) given 

in terms of, either a sequence of Wiener processes (see, e.g., Einmahl, 1992 and 
Dress et ah, 2006) or a single Wiener process as in Proposition 3.1 of Einmahl et al. 
(2006), will be very crucial to our proof procedure. In the sequel, we use the latter 
representation which says that: there exists a Wiener process W, such that for every 
0 < 7 < 1 / 2 , 

sup \an (s) — W (s)| A 0, as u —)■ cxD. (5-12) 

0<s<l 

We begin by hxing xq > 0, then we decompose k~^^‘^'Dn (x ), for x > xq, as the sum 
of the following four terms: 

_l/7l (^xXn—k:n) F {xX^—k-.n^ 

F {xXn-k-.n) 


M„1 (x) := X 
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P = 

0.7 







7 i = 0.6 




o 

bo 


N 

n 

k* 

absolute bias 

rmse 

n 

k* 

absolute bias 

rmse 

200 

139 

8 

0.1811 

0.4645 

140 

9 

0.2485 

0.6034 

300 

210 

17 

0.1280 

0.3451 

209 

23 

0.1230 

0.5082 

500 

348 

28 

0.1151 

0.2803 

348 

31 

0.1024 

0.3708 

1000 

699 

43 

0.0461 

0.2421 

700 

42 

0.0684 

0.3825 

1500 

1050 

67 

0.0212 

0.2362 

1050 

57 

0.0527 

0.2539 

2000 

1399 

95 

0.0254 

0.2261 

1401 

94 

0.0469 

0.2602 





P = 

0.8 




200 

159 

10 

0.0759 

0.5235 

159 

11 

0.1435 

0.6815 

300 

240 

21 

0.0485 

0.3647 

238 

18 

0.1052 

0.4511 

500 

399 

34 

0.0403 

0.2718 

400 

41 

0.0935 

0.3170 

1000 

800 

55 

0.0408 

0.1999 

800 

54 

0.0504 

0.2881 

1500 

1205 

116 

0.0562 

0.1911 

1199 

85 

0.0428 

0.2362 

2000 

1599 

117 

0.0285 

0.1534 

1599 

117 

0.0373 

0.2054 





P = 

0.9 




200 

180 

14 

0.0537 

0.5204 

179 

17 

0.1098 

0.7531 

300 

269 

24 

0.0388 

0.3343 

268 

26 

0.0844 

0.3654 

500 

450 

48 

0.0294 

0.2869 

450 

49 

0.0721 

0.2448 

1000 

899 

64 

0.0359 

0.1557 

899 

67 

0.0465 

0.2014 

1500 

1349 

103 

0.0171 

0.1350 

1349 

101 

0.0385 

0.1828 

2000 

1799 

144 

0.0188 

0.1107 

1799 

145 

0.0251 

0.1466 


Table 4.1. Absolute bias and raise of the tail index estimator based 
on 1000 right-truncated samples of Burr models. 


M „2 {x) 


F {xXn—k:ri) Fn, F (^Xn—k:n} 

Fn (^n—k:n) F (^n—k:n) 


^ F ^xXfi—k:n) Ffi (xXn—k:n) F ^xX^—k-.n) 

\ Fn (yX^—k'.n) ) F {xXn—k'.n) 


and 


M„4 (x) 


F {xX^i—k-.n) 

^- X ' 

F {Xn-k:n) 
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In order to establish the result of the theorem, we will successively show that, under 
the hrst-order of regular variation conditions, we have uniformly on x > Xq, for 
7/72 <17 < 1/2 and e > 0 sufficiently small 

(x) 

= x^!^ I^W [x-^!^) + ^ ^ t-7/72-iw dtj + Op (e) 

\/fcM„2 {x) = -[ — W ( 1 ) + w dt\+Oj, (e) 

17i 7i + 72 Jo j 

and 

(x) = Op (e) x-V7i+(i-^)/7±y 

Moreover, if we assume the second-order condition we will show that 

"^1 /'yi T 

\/fcM „4 (x) = (1 -I- Op (1)) - - VkAo in/k). 

llTl 

Here Op and Op stand for the usual stochastic order symbols. For convenience, let 
Ofc := \]p{n/k) and recall that Ui? is regularly varying (with index 7 ). Then by 
combining Corollary 2.2.2 with Proposition B.1.10 in de Haan and Ferreira (2006), 

p 

we show that Xn-k-.n/cbk —)■ 1 as n —)■ 00 , which implies, due to the regular variation 
of F, that F (xttk) /F (xXn-k-.n) = 1 + Op (1) and therefore 

M„i(a;) = (l + Op(l))M;i(a;), (5.13) 


where 


]y[* _ jxXn-km) F (xXn-k-.n) 

F{xak) 


Now, observe that, in view of equation (1.5), we may write 


F (x) = exp —A (x) and F„ (x) = exp —An (x ), 


where A (x) and its empirical counterpart A„ (x) are dehned by dF (z) jC {z) and 
/“ dFn (z) jCn {z) respectively. Note that F„ {xXn-k-.n ), F {xXn-k-.n) and F (xa^,) 
tend to zero in probability, uniformly on x > xq, it follows that A„ (xXn-k-.n), 
A {xXn-k-.n) and A (xa*,) go to zero in probability as well. Using the approximation 
1 — exp(—f) ^ t, as 11 0, we may write 


{ x ) = (1 + Op ( 1 )) 


An (xXn—k-.n) A {xXn—k-.n) 


A (xofc) 

Next, we provide a Gaussian approximation to the expression 


/~rAn {xXn—k'.n) A {xXn—k-.n) 

Aj^) 
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then we deduce one to (x). For this, we decompose the difference 

{xXn-k-.n) - A {xXn-k-,n) ^ the sum of 

r^^Fn{z)-F{z)) 


Snl(x) := - 


C{z) 


S „-2 (x) := - 




dFfi , 


and 


Sn3 (x) := 


rxak 


Cn (z) C (z) 
d{F^{z)-F{z)) 


JxX„-k:n ^ 

For the hrst term, we use the fact that Fn (z) = 0 for 2; > Xn-.n, to write, after an 
integration by parts and a change of variables, Sni (x) = (x) — (x), with 

(X) := T'n(a.x) -F {a,x) f „ (g^) - F (a,.) ^ 


di’ (x) : = 


C (Ofcx) 


C2 (ttkZ) 


y/k /n— 


It is easy to verify that Fn (xak) — F (xak) = -( —F (xak) ] , it follows that 

n \k 


y/kSni {x) w 

= dn (x) an yF (OfeX) ) , 


A (Ofcx) 


where dn (x) := 


k/n 


From Lemma 6.2 {in ), we have 


A (xofc) C (ofcx) 

dn (x) = (7/71) x^/'^ + O (e) 


(5.14) 


as n —>■ oc, uniformly on x > Xq, it follows that 
x/fcAlV (x 


= {( 7 / 71 ) x^'^^ + Op (e)x^/'^=^"} {^F{akx) 


A (ofex) 

On the other hand, for 0 < 7 < 1/2, the sequence of rv’s supo<g<]^ |a„ (s)| /s^ is 
stochastically bounded. This comes from the inequality 

sup \an (s)| < sup S~'^ \an (s) — W (s)| + sup s~'^ |W (s)|, 

0 <S<1 0<S<1 0<S<1 

with approximation (5.12) and the fact supo<5<i s”’'|W (s)| = Op (1) (see, e.g.. 

Lemma 3.2 in Einmahl et ah, 2006). Now, let e > 0 be sufficiently small. Then, 

by applying Potter’s inequalities to F (see, e.g.. Proposition B.1.9, assertion 5 in 

de Haan and Ferreira, 2006), we write —F{akx) < (1 + e)x“^0±e^ it follows that 

k 

an {(^kx)^ = Op _ por notational simplicity and without loss of gen¬ 

erality, we attribute e to any constant times e and to any linear combinations of 
y±cie i;=‘='=2^^ for every Ci, C2 > 0. Therefore 

VkSS jx) _ 7 ^1/. , fn- 


A (ofcx) 


= -X ' 'ttr 

7i 


k 


F(afcx)) +Op (6)x(i-")/^±^ 
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For (x), let us write 

C{aux) /■“ C\au) (n 
^{aux) C{ak) j, C^{auzf^\k 


C{akz) 
C{ak) ■ 


From Lemma 6.2 (i ), the function C is regularly varying at infinity with index 
(—1/72) , as G is, this implies that C {xak) jC {ak) = + O (e) Then 

by using (5.14) , we get 


dn (a:) 


C (ttkx) 

C{ak) 


(7/71) + O (e) as n —>■ 00. 


(5.15) 


For convenience, we set (x) /A (a^x) = {1 + O (e) x^*^} Tn (x), where 

6*2 (afc) 


Tn (x) := 

7i 


(Xr 


(ttkz) \k 




which we decompose in the sum of 


In (^) • 


- 


(ak) 


7i 


(ttkz) \k 


n 


an -rF (ttkz)) d 


C{ak) ’ 

F jakz) 
Ciak) ’ 


Jn (2^) 


-XxFzi 

7i 



_ ^2/72 

(akz) 


CXr. 


77_ 

-F {ttkz) ] d 


G (ttkz) 
G{ak) ’ 


and 

Kn (x) := —x^Fi f (^p (^akz)) dyypy. 

7i Jx / G(ak) 

Recall that Uk oc, G (ak) ~ G (ak) and F (ak) = o (G (ak)) as n —)■ 00. On the 
other hand, by using, once again. Potter’s inequalities to O, (regularly varying at 
inhnity with index —1/72), we write, for all large n and z > x, 


(1 — e) z (p z ’^) < 


G (akz) 
G(ak) 


< (1 + e) z-F'^^max {z%z-^) . 


(5.16) 


It is clear this implies that (ak) jG’^ (akz) < (1 — e) ^ g|.Q_ 

chastic boundedness of supo<s<i |a„ (s) | /s’’ and the fact that —F (akz) < (1 + e) 

~ fC 

we have 

dn (2:) = Op ( 1 ) X^F^ f ^2/72^2. (z-Fz±ey ^ d^i^kZ) ^ 

^ ’ J. ^ ’ F(ak) 

Integrating by parts, we readily get T (x) = Op (1) a;V72-r?/7±e = ^(i-??)/ 7±7 pet 

us now consider Jn (x). From Proposition B.1.10 in de Haan and Ferreira (2006), we 

have \G (akz) /G (ak) — z~^F 2 1 < ez ^h 2 ±<^_ Applying the mean value theorem, then 

combining this inequality with (5.16) , yield 


(^k) _ 2/72 <- 2 (z^'' + 1 ) 2h2±e 

G^akz) - (l-ef 
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Similar arguments as the above lead to (x) = Op (e) »?)/7±e_ Now, we focus on 

Kn {x). Since C {ak) ^ G {ak), then 

(x) = (1 + Op (1)) —(yF (akz)) 


7i Jx } G[ak) 

Let denote the quantile function pertaining to df G and use the change of 
variables (l — sG (ak)) /a^ to get 


G(afcai) / ^ \ 2 Ho 

^ 1/71 f GK) f G"^ (l — sG (qfc)) A 


K^{x) = -{l + Op{l))^x^Fi / 

7i Jo 


_ _ 

where in (5) := —F (l — sG (a^))) . It is easy to check that 


*^n (^n (^)) 


K^{x) = -{l + Op{l))Y,K^ 


where 


Gjakx) 


ix) := ^x^ [ f 


7i Jo 


^ (yin (^)) ds^ 


G(akx) 


Kn 2 {x) := —X^Fl / _5 ^an{in{s))ds, 


Knsix) := —x^Fi / s ‘^an{in{s))ds. 

7i Jo 

By routine manipulations and similar arguments based on stochastic boundedness of 
supo<5<i Icin (<s)| /s^ and the aforementioned Proposition B.1.10 applied to the regu¬ 
larly varying functions G and G^ (1 — •), we show that Kni (x) = Op (e) 
i = 1,2 and iL„3 (x) = Op (1) therefore we omit the details. Up to this 

stage, we have shown that Tn (x) = Op (1) It follows that 

which, after gathering the components of Tn {x ), is equal to 


—x^Fi / 

7i Jo 


s ^an{in{s))ds + Op{e)x^^ 


Therefore 


'^Snl (X) _ 1 „ C^p, 

A (an) T, Q„qj-(an)j 


— —x^Fi / 

7i Jo 


S “^an {in {s)) ds + Op {e) X^^ »?)/7=Fy 
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Recall that 71 < 72 and 7/72 = 71/ (71 + 72), then we may choose the constant 7 
in snch a way that 7/72 < 7 < 1 / 2 . Making nse of weak approximation ( 5 . 12 ), we 
obtain 

VkSni jx) 

A {akx) 

»3;-1/72 

= (ofex)) + ^x^/^ / s-^W (4 (s)) ds + Op (e) xA-^^^O^h 

7i / 7i io 

Note that k/n = F (O^ (l — G (0^))) , hence 

F(G^(l-sG(a,))) 

F(G-(1-G(a,))) ■ 

Since s —>■ FoG^ (1 — s) is regularly varying at infinity with index 72/7, then, from 
Proposition B. 1.10 in de Haan and Ferreira ( 2006 ), we have for all large n 

oon (s) := \in (s) — ( 5 . 17 ) 


Recall that Xq > 0 is fixed, then sup2,>,j,^ suPq^^< 3 ,-i /72 oJn (s) —)■ 0, as n ^ cx). On 
the other hand, by using Levy’s modulus of continuity of the Wiener process (see, 
e.g.. Theorem 1.1.1 in Csorgo and Revesz, 1981), we have 


W (4 (s)) - w (4^/^) I < 2^0Jn (4 log(l/a;„ (s)). 


uniformly on s > x“^A2^ almost surely. By using the fact that, logu < as u 0, 
together with inequality (5.17), we get |W (4 (s)) — W | < 

Following our convention, we may write that ( 72/7 ± e) (1 — e/2) = 72/7 ± e. Since 
7i < 72 then 72 / (27) > 1 and after elementary calculation, we show that uniformly 
on X > xo 


7i 


C-l/72 


— A ^1/71 


C-1/72 


s (4 (s)) ds = —x^ 
7i 


s ds+Op (e) x^A 27)±7 


By similar arguments, we get 


—x^Aw (yF {akx)) = ^x^/^W (x"^/^) + Op (e) x^A 27)±7 
7i ^k / 7i 

It is obvious that Op (e) x^A27)±e _|_ Op (e) xA“^)/'i'=‘=^ = Op (e) xA“^)4=‘=^, it follows 
that 


VkSni jx) 
A (ofcx) 


= —x^Aw (x-^/^) + 
7i 


—x^/^i 
7i 



s-2w (4^/^) ds + Op (e) xA-^)4±7 






15 


After a change of variables, this may be rewritten into 
\/kSni (x) 


A (akx) 


(5.18) 




= —X- +- [ — 3 . 1/7 / ^-7/72-1^ (x-^AA fit + n (e) 

71 ^ ^ 71 + 72 io ^ ’ 

Now, we consider the second term Sn 2 ix). We have {z) = 0, for 2; > X„.„, thns 


rXr, 


Therefore 


Sn2 (X) = 


\Sn2 (ai)! < 0, 


Cn{z)-C{z) 


C^{z)C{z) 

r mz)-c{z)\ 


dFn (z). 


dFn (z ), 


JxX , C'2 (z) 

where 9^ := snpj^^ {C (^) /C'„(^)}, which is stochastically bonnded (see, 

e.g., Stnte and Wang, 2008). By recalling that C = G — F and Cn = Gn — Fn, with 
Gn denoting the empirical df of G, we write |S'„2 (a^)| < dn (T„i (x) + T„2 (a:)), where 

|F„(z)-F(z)| 


and 


Tni (a:) := 


Tn 2 (a:) : = 


Lx , L) 

p \Gn{z)-G{z)\ 
G2 iz) 


-dFn (z) 


■dFn {z). 


— k:n 

It is easy to verify that, by a change of variables, we have 
VkTni (x) , / N k/n G (ofcx) 


A (ofcx) 


= dn (a:) 


G{ak) G{ak) 


X 


G2 (afc) 


^ F [xXn-k-.nZ]) {xXn-k:nZ) 


G 2 (xX„_fc:„) G 2 (xXn-k-.nZ) jG‘^ {xXn-k,r^ ‘ F (Ofe) 

Recall that, nniformly on x > xq, we have G (ak) IG (xXn-k-.n) = Op (1) 
Moreover, we nse (5.16) and (5.15) to write 

k/n 


VkTni (x) _ 

A(afcx) ^\G{ak) 


X 


l/ 7 ±e 


X 


^ 2/72 


Cln {xXn—k:nZ'^ 


d 


Fn {xXn—k:nZ^ 

F {ak) 


On the other hand, by nsing the stochastic bonndedness of snpo<5<i \an (s)| /s^ we 
get 


y/kTnl (x) ^ ^ f k/n \ f 2 h 2 -rth±eA ^ri (xXn-h.nZ) 

A{akx) ^^\G{ak)) A F(a,) 
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where the integral may be split as follows 


F{ak) 


where 


and 


Pn (x) := / 


2/'y2—??/7±e^ j Fn ipzXn—k-.nF) F (^xX^—k-.n^') 


F{ak) 


Qn (x) := j z 


2 / 72 -» 7 / 7 ±« 7 -^ i^Xn-k-.nZ) 

F{ak) 


It is clear that 

(x) = k-F^ (^F (xX„_,.„z)) . 

By similar argnments as those nsed above, we show that F„ (x) = Op 
Qn (x) = Op {x~^F^<^'^ _ Therefore 

^Tnl (X) ^ 

A(afcx) 

Next, let Vi \= G ( 1 ^), z = 1,..., n, and define the corresponding tail empirical process 
fdn (s) := Vk (V„ (s) - s) , for 0 < s < 1, where V„ (s) := Yh=i 1 (Fi < ks/n) . 
Like for a„ (s), we also have snpo<s<i |/5n (s)|/s^ = Op (1), therefore by similar 
argnments as those nsed for T„i (x) , with the facts that G (t) ~ O (t) as t —)■ cxd and 
72 > 7, we show that 

(X) f k/n \ 

A(a,x) ^\C^-^iak)) 

Tt Tl 

From Lemma 6.2 (ii), we have that both —G {ak) and {ak) tend to infinity, 

fC fC 

it follows that 

VkTnl (x) 


'^'F'n2 (x) ^ 

A(afcx) ^ A(afcx) ^ 


= Op (x and 


x^ 


) 


Since Op (x ^0='=^) + Op (x^ ’?)/7±ej = (a;(i f;)/7±ej ^ then 


\/kSn 2 (x) 


= Op lx 


(l-»7)/7±e 


) 


(5.19) 


A (ofex) 

Let ns now focus on the third term F„3, which, by integration by parts, equals the 
sum of 


Snl {x):=- 


f*xa)^ 


— k.-.'n 


Fn {z)-F{z) 
G2 G) 


dC {z ), 


(2) _ Fn (OfcX) F {cikX') ^ Fn {xXn—k-.Q F (^xXn—kin) 


(x) = 


G (ofcx) 


+ 


G {xXn—k'.n) 


and 
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By using the change of variables 2; = txa^ we get 


(x) 


A {a^x) 


= -dn (x) 


an[^F{akXz)^ C {a^xz) 


and 


^/7sS' 


lx 


A (akx) 


= -dn (x) <{ ( -F (ofcx) 


n 


Xn-k:Jak (C (akxz) /C{akx)f C (a^x) ’ 
C (ofcx) 


C{xXn-k:n) Vfc 


n 


( > F ip^^n—k'.n) 


Routine manipulations, including Proposition B.1.10 in de Haan and Ferreira (2006) 
and the stochastic boundedness of supo<s<i (s)| /s^, yield 




[X) 


A (ofcx) 


A (ofcx) 

It follows that 

(x) /A (Ofex) = Op 

By gathering results (5.18) , (5.19) and (5.20) , we obtain 
pr^n {p^^n—k'.n^ A (^xXn—k:n) 

^ AK^) 


= Op (X 


(l-r?)/7±e 


)• 


(5.20) 


(5.21) 


= (x-^/^) + -XL — (x-^FA dt + Op (e) 

7i ^ ^ 7i + 72 io ^ ’ 

which yields that 


(x) 

= xF'^ I^W {x-F^) + ^ ^ + Op (g) a;(i-^)A±7 

We show that the expectation of the absolute value of the first term in the right- 
hand side of the previous equation equals Op (x^A27)j _ Since 1/ (27) < (1 — p) /y, 
we have x^Aiy^M*;^ (x) = Op , which leads to 

\/fcM„i (x) = \/fcM* 1 (x) -I- Op ^ 


Recall that e > 0 is chosen sufficiently small, then for any 0 < p < 1/2, we have 
\/fcM„i (x) 

= x^A {x-F'^) + ^ ^ ^-7/72-1^ (x-^A^) + Op (e) xA-^)A±7 

Before we treat the term M„2 (x), it is worth mentioning that by letting x = 1 in 
the previous approximation, we infer that 

_ 1 = Op {k-p = Op (1). (fc ^ DO) . (5.22) 

/c:nj 














18 


This, with the regular variation of F, imply that 




= (l + Op X 




Fn {Xfi—k-.n} 

To represent \/fcM„2 (ai), we apply results (5.21) (for x = 1) and (5.23) to get 


(5.23) 


(a:) = - <1 — W (1) + 

.71 71 + 72 


7 


^-7/72-1^ (f) dtl + Op (e) 


X 


For the third term M„3 (x), we write 


(a:) = x-^/A x'/^l^/fcM„l (x), 

which, by equation (5.23) , is equal to Op (e) a;“^Oi+d-+/7±q pet r/o be such that 
7/72 < ho < 7 < 1/2, then rjo — rj < 0 and for e > 0 sufficiently small, we have 
(ho — h) /7 + e < 0. Since x > Xq > 0, then Op (e) = Op (e) and thus 


3.i/7i-(i-r,o)/7 (M^i (x) + M„2 (a;) + M„3 (a:)) - F (x; W)} = Op (e), (5.24) 


where F (x; W) is the Gaussian process given in Theorem 2.1. For the fourth term 
M„4 (x) , it suffices to use the uniform inequality to second-order condition (2.9) , 
given in assertion (2.3.23) of Theorem 2.3.9 in de Haan and Ferreira (2006), to get 

VkMnA (x) = (1 + Op (l))x~^'^'^^- - -\/kAF (l/F (Xn-k:n)) , 

7i'^i 

for a possibly different function Ap with Ap ~ Ap. Then Proposition B.1.10 in 
de Haan and Ferreira (2006) and the fact that t -7 Ap (l/F (Up (t))) =: Aq (t) is 
regularly varying with index ri/71 with Xn-k-.n/cLk ^ 1, imply that 


Aq {l/F (Xn-k-.n)) p 

-T—=-r->■ 1, as n —)■ CX7, 

A„(1/FK)) 

2;Fi/71 _ I 

as well. Since Op(l)x“^+i- = Op , and by assumption 

VkAo [l/F {ttk)) = VkAo {n/k) = O (1), it follows that 


^ , q-ri/jl _ 1 

VklSAni (x) = x~^Fi - 


7iU 


\/fcAo {n/k) + Op (x _ 


Finally, by letting e 0 in (5.24) , we end up with 


sup X 

X>Xq 


l/7l-(l-»7o)/7 


E)n (x) — F (x; W) — X 


— -r- ^+1 


Xn /71 _ 1 


7iU 


VkAo {n/k) ^ 0, 


for every xq > 0 and 7/72 <ho<h<l/ 2- Letting 170 := 1/2 — and recalling that 
1/71 = 1/7 — 1/72 yields that 0<^< 1/2 — 7/72 and achieves the proof. □ 
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5.2. Proof of Theorem 3.1. We start by proving the consistency of 71 that we 
write as 71 = x~^Fn (xXn-k-.n) /Fn (Xn-k-.n) dx. It is readily checked that this 
may be decomposed into the sum of 

3 


— 1 ^ {p^^n—k'.r 


Iln := / X 

A F {Xn-k-,n) 


-dx and 


/ OO 

x~^ ^ 

i=i 


(x) dx. 


By the regular variation of F (1.1) and Potter’s inequalities, we get —)■ 71 as 

n —)■ OO. Then, we just need to show that l 2 n tends to zero in probability. From 
(5.24) we have 

1 POO -1 POO 

J2„ = — / x-^T (x; W)dx+^ / x-^Op dx. 

yk Ji yk Ji 

On the one hand, since 7/72 < 7 , the second integral above is finite and therefore 
the second term of l 2 n is negligible in probability. On the other hand, we have 


POO POO 

/ x-^F (x; W)dx=— / {W (x-^/^) - x-^/^W (1)} dx + 

Ji 7i Ji 


7 


7i + 72 


POO 


1/72-1 J / g-7/72-1 {s)} ds \ dx. 


which, after some elementary but tedious manipulations of integral calculus (change 
of variables and integration by parts), becomes 


/ OO 

x-^r(x;W)dx = -7W(1) 


(5.25) 


7 


/ (72 — 71 — 7 log s) s (s) ds. 

7i + 72 Jo 

By using the facts that E | W (s) | < and 71 < 72 , we deduce that x“^r (x; W) dx 
is stochastically bounded and therefore the hrst term of l 2 n is is negligible in prob¬ 
ability as well. Consequently, we have l 2 n = op (1) when n ^ 00 , as sought. 

As for the Gaussian representation result, it is easy to verify that y/k (71 — 71 ) = 
/j°°x“^D„ (x) dx. Then, applying Theorem 2.1 yields that 

^ + r x-^T (x; W) dx + op (1), 

1 - !■ Jl 

and finally, using result (5.25) completes the proof. □ 


5.3. Proof of Corollary 3.1. We set 


(71 - 7i) = 7^ + + Op (1), 


1 — r 
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where A := aAi + 6A2 - A3, with a := (72 - 71) / (71 + 72) , & := -7/ (71 + I2) and 

Ai := f sP-^W{s)ds, A2 : = 

Jo 

with p := 1 — 7/72 > 0. 

It is clear that the asymptotic mean is equal to limjj_j.oo V^Aq (n/k) / (1 — r), while 
for the asymptotic variance we hnd, after elementary but tedious computations, the 
following covariances: 


f (s) log sds, As := W (1), 

Jo 


E[A?] 
E [A1A2] 


2 

P(2p-1)’ 


E[A^] 


2(4p-l) 

{2p-ir 


E[A^] 


1 , 


1 -4p 
(2p-l)^’ 


E [Ai As] 


1 

P 


E [A2A3] 


1 

p2' 


It follows that 


F rA21 = 2a^ 2b^ (4p - 1) 2a&(l-4p) , ^ ^ . 

Replacing a, b and p by their values achieves the proof. 


□ 


Concluding notes 

We would like to emphasize the fact that, unlike Gardes and Stupfler (2015) who 
defined their estimator in terms of two (not necessarily equal) sample fractions 
k = k' oi upper order statistics from X and Y respectively, we build our estimator 
on the basis of just a single sample fraction. The consideration of two distinct sample 
fractions poses a problem from a computational point of view, as the issue of selecting 
an optimal couple of sample fractions is not as easy and usual as determining just one 
best number of top statistics to be used in parameter estimate computation. Besides 
that. Gardes and Stupfler (2015) didn’t treat the asymptotic normality when k = k' 
and only carried out their simulations in this very particular case, as they mentioned 
in their conclusion. For these reasons, we don’t compare the two estimators in 
Section 5. 


A more thorough simulation study, with conhdence interval construction and even¬ 
tual comparison with the estimator of Gardes and Stupfler (2015), will be part of a 
future work. Another point, beyond the scope of the present paper, that deserves to 
be considered is to reduce estimation biases under random truncation. Similar an¬ 
terior works were done with complete datasets by, for instance, Peng and Qi (2004), 
Li and Peng (2010) and Brahimi et al. (2013). 
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We finish this work by making a comment on relation (3.11) , which actually is a 
special case of a more general functional of the distribution tail dehned by 



where g is some weight function and a some positive real number. As a consequence 


of the fact that limt^oo h* {g, a) = 7 ", this functional can be considered as the start¬ 


ing point to constructing a whole class of estimators for distribution tail parameters. 
Indeed, in the complete data case, we replace F by its empirical counterpart and 
t by Xn-k-.n to get the following statistic which generalizes several extreme value 
theory based procedures of estimation already existing in the literature: 



When = a = 1, we recover the famous Hill estimator (Hill, 1975). For a detailed 
list of extreme value index estimators drawn from the statistic above, we refer to the 
paper of Ciuperca and Mercadier (2010), where the authors propose an estimation 


approach of the second-order parameter by considering differences and quotients 


of several forms of Tn,k {g, «) • By analogy, when we deal with randomly truncated 
observations, we substitute the product-limit estimator F„ for F in the formula 
of rt{g,a) in order to obtain the following family of parameter estimators under 
random truncation: 



This would have fruitful consequences on the statistical analysis of extremes under 
random truncation. 
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6 . Appendix 

Lemma 6.1. Assume that both second-order conditions (2.9) and (2.10) hold. Then, 
for all large x, there exist constants di,d 2 > 0 , such that 

F (x) = (1 + o ( 1 )) dix~^^'^ and G (x) = (1 + o ( 1 )) d 2 X~^^'^^. 


Proof. We only show the hrst statement since the second one follows by similar 
arguments. To this end, we rewrite the hrst equation of (1.2) into 


F (x) = —p (x) F (x 


G {xz) {xz) 
L G(x) F(x) 


By applying Proposition B.1.10 in de Haan and Ferreira (2006) to both F and G, 
it is easy to check that 


G {xz) {xz) 
1 G (x) ^ F (x) 


(l + o(l))7/7i. 


On the other hand, since F and G satisfy the aforementioned second-order condi¬ 
tions, then in view of Lemma 3 in Hua and Joe (2011), there exist two constants 
01,02 > 0, such that F (x) = (1-|-o (1)) OiX“^/'^i and G (x) = (1-|-o (1)) 02 X“^/'^ 2 , 
as X ^ cx). Therefore F (x) = {1 o (1)) dix~^^"' with di = ^“^ 01027 / 71 . □ 


Lemma 6.2. Under the assumptions of Lemma 6.1, we have 
(i) lim C {t) /G{t) = 1. 

t—>-oo 

(a) lim (Uf (t)) = 00, for each 0 < u < 1. 

t—)-CxD 

{Hi) limi^oo sup3,>^p \ {tA (xU^ (t)) C (xU^ {t)))~^ - (7/71) x^G\^ = 

for xo > 0 and any sufficiently small e > 0. 

Proof. For assertion (i), write C {t) /G {t) = 1 — F {t) /G {t) and observe that from 
Lemma 6.1 we have F (t) /G (t) = (1 -|- o (1)) {di/d 2 ) pF 2 -^h^ Since I /72 — I /7 < 0, 
then F {t) /G{f) = o(l), that is G{t) /G{f) = 1 -|- o(l) as sought. For result 
{ii ), Lemma 6.1 implies that Up {t) = (1 -l- o (1)) {dpy (as t —)■ cxd), it follows that 
G (Ui? (t)) = (1 -I- o (1)) d 2 {dit)~'^^'^^. Since 0 < 7/72 < 1, then for every 0 < i/ < 1, 
P^G (Ui? (t)) —)■ cxD as f —)■ 00. To prove {Hi ), we hrst show that 

tA (xUi. {t)) G (xUi. (t)) - (71/7) x-^/^ = O (e) (6.26) 

Recalling that A (x) = dF {z) /G {z) and F {Up {t)) = 1/t, we write 

G(xUi.(t)) G(Ui.(t)) dF{zl]p{t)) 

C'(Uf {t)) y, G{zl]p{t)) F{iJp{t)) 


tA {xVp {t))G{xVp {t)) 


(6.27) 
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Observe now that tA (xUi? (t)) C (xUi? (t)) —may be decomposed into the 

7 

snm of 


f CjxVF (t)) r CjlJFit)) dFjzVFit)) 

[c{iJF{t)) "" ) C{zVF{t)) F(U^(t)) 


D2 (s;t) 





( C{VF{t)) _ 
(t)) 


dF{z\]F (t)) 
F{VF{t)) 


L>3 (s;t) : = 


•= —X~^F'2 


and _ 

V f (UfW) 

By applying Proposition B. 1.10 in de Haan and Ferreira ( 2006 ) to both C and F 
with integrations by parts, it is easy to verify that 


\tA (xUf (f)) C (xUf (t)) - (71/7) x-F'^\ < 

Observe now that tA (xUi? (t)) C (xUi? (t)) — (7/71) x^F jg eqnal to 
((fA(xUF(t))C'(a;UF(t)))"^) ^-((71/7)2:"^/^) ^ 


By nsing the mean valne theorem, the latter eqnals 

(71/7) X~^F — (xUiT (t)) C (xUiT {t)) 

where xf) {x; t) is between (71/7) x~^F am^ tA (xUi;’ (t)) C (ajUF {t)) ■ In view of the 
representation ( 6 . 27 ) and Potter’s ineqnalities, applied to C and F, with an inte¬ 
gration by parts, we get tA (xU_p (t)) C (a;U_p (f)) > (1 — e) it follows that 

{;ijj {x]t))^ > (1 “ e)^ and therefore 

|fA (a:U_p [t)) C (xU_p {t)) — (7/71) x^F^ < (1 — e)~^ 


as songht. 


□ 








